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We study the effective potential of a real scalar ip 4 theory as a function of the temperature 
T within the simplest ^-derivable approximation, namely the Hartree approximation. We apply 
renormalization at a "high" temperature T* where the theory is required to be in its symmetric 
phase and study how the effective potential evolves as the temperature is lowered down to T = 0. In 
particular, we prove analytically that no second order phase transition can occur in this particular 
approximation of the theory, in agreement with earlier studies based on the numerical evaluation 
or the high temperature expansion of the effective potential. This work is also an opportunity to 
illustrate certain issues on the renormalization of 5>-derivable approximations at finite temperature 
and non-vanishing field expectation value, and to introduce new computational techniques which 
might also prove useful when dealing with higher order approximations. 
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I. INTRODUCTION 

The two-particle-irreducible (2PI) effective action [IHE] offers a general and systematically improvable approach 
\ for resumming infinite classes of Feynman diagrams of a given quantum field theory. One of its most compelling 
©H 1 aspects is that it appears to be applicable to a wide variety of situations ranging from out-of-cquilibrium settings to 
(-^ | finite temperature calculations in equilibrium. In fact, the first orders of approximation of the 2PI effective action 
. seem to capture already many interesting features of quantum field theories. In the case of scalar fields for instance, 
' including the first non-local, field independent contribution to the 2PI effective action is enough to obtain a controlled 
(non-secular) time evolution which shows thermalization at late times, at least at the level of the two-point function 
[6(. In equilibrium, the same orders of approximation lead to a rather good convergence of some thermodynamic 
| quantities, such as the pressure [3], indicating that part of the infrared physics is properly taken into account. These 
Q\ ' observations are not limited to scalar theories but concern also theories coupling scalar and fermionic degrees of 
freedom (out-of-equilibrium) Q or gauge theories (at finite temperature) 0,^3] . Some more genuine non-perturbative 
effects can be captured by combining the 2PI effective action with a large-N type expansion [111, [l2j ■ In the non- 
equilibrium context the 2PI 1/N expansion at next-to-leading order proved to be particularly fruitful for the study of 
C*") , a variety of problems, such that thermalization fliT - [l3T ] , preheating [lj, [Till , transport coefficients [l6l - fl8| , non-thermal 
fixed points [l!|, decoherence [2(| and topological defect formatio n [2 ill . In equilibrium, the 2PI formalism was 
applied to phenomenological studies in various approximations, see [22| and references therein as well as [23|,[24[. Of 
course, one can easily imagine that certain non-perturbative features are beyond reach within the 2PI framework or 
its present approximations. However, hig her order approximations of the 2PI effective action or its generalization 
to 3PI, 4PI, NPI effective actions [25(, although numerically unaffordable at present, could allow to capture an 
increasing number of such effects in the future (26l - |28j . Some investigations in these directions can be found in (29j - [3l| . 
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Concomitantly with the increasing number of applications of the 2PI effective action, new insight has been gained 
on technical aspects regarding its renormalization. The systematic renormalization of approximations based on the 
2PI effective action was first understood in the case of scalar theories at finite temperature both in the real-time 
formalism [32l . [33| and in the imaginary time formalism |3~H. l35l . It was then extended to the case of scalar theories 
coupled to fermions (36| and also to abelian gauge theories |37l |. Another important extension consisted in including 
the possibility of non-vanishing field expectation values, of relevance for the determination of the effective potential 
and thus the study of phase transitions (3f| I38l - l42l ] . The renormalization of models with more complicated global 
symmetry was studied in (43[. From all these studies, a certain consensus has emerged that the renormalization 
of approximations based on the 2PI effective action seems always possible provided certain extensions of the 
renormalization procedure are allowed in order to cope with certain approximation artifacts. For instance, in the 
presence of non-vanishing field expectation values, it is well known that there exist different expressions for the two- 
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and four-point functions, and more generally for higher n-point functions. Although equivalent in the exact theory, 
these various expressions differ within a given approximation and bring their own divergences which can only be 
absorbed by allowing for apparently more counterterms than usual [35|. Some of these counterterms are fixed by 
means of the usual renormalization conditions. The others are fixed by imposing consistency conditions, that is 
conditions which would be satisfied automatically if no approximation was considered at all [3a . [iij . The fact that 
the consistency conditions do not involve the parameters of the theory is crucial to maintain the number of such 
parameters to its expected value, despite the larger number of counterterms. In fact, these apparently different 
counterterms should be viewed as different subseries of the complete perturbative countcrterm series and thus should 
agree when no truncation of the 2PI effective action is considered. 

In this paper we illustrate partly these issues regarding renormalization by revisiting one of the simplest approx- 
imations of the 2PI effective action at non-zero field and finite temperature in the case of a real scalar field, the 
so-called Hartree approximation. The calculation of the finite temperature effective potential of a real scalar ip A 
theory in the Hartree approximation and the order of the phase transition that it predicts have been discussed many 
times in the literature O EH-EH- Similar discussions exist for different approximations in the O(N) model with 
N > 2, see [4!|. There is a wide consensus on the fact that the order of the phase transition predicted by this 
approximation is first order, even though all existing studies are either based on numerical evaluations for particular 
values of the parameters or involve a high temperature expansion which is not justified for all values of parameters. 
Part of the originality of this work is that, without relying on a high temperature expansion, it presents a complete 
analytic confirmation of these results, in the whole parameter space. 

In section [Til wc briefly recall how the effective potential is computed using the 2PI effective action. The main 
ingredient is the resolution of a self-consistent equation for the two-point function, the so-called gap equation, which 
encodes the resummation of particular classes of Feynman diagrams. The Hartree approximation corresponds to the 
leading order of a systematic approximation scheme for evaluating the 2PI effective action and the corresponding gap 
equation is then an equation for a self-consistent mass. This equation encodes the resummation of superdaisy diagrams 
\m , IHol . l5lj and is discussed at length in section IIIII We first recall how renormalization of the self-consistent mass is 
performed at finite temperature and non- vanishing field expectation value, following a slightly different presentation 
than the one which is usually found in the literature (although equivalent in practice). We put the emphasis on the 
fact that not all the divergences of the gap equation lead to divergences of the self-consistent mass. The quadratic 
divergence of the gap equation is indeed responsible for a divergence of the self-consistent mass, which needs to be 
absorbed in the renormalization of the bare mass, as usual (sections MI Al and MI B[) . In contrast, the remaining 
logarithmic divergence of the gap equation does not lead to a divergence of the self-consistent mass. Instead, it makes 
this self-consistent mass trivial in the continuum limit, that is equal to the renormalized mass, independently of the 
value of the temperature or the field expectation value. In this respect, the renormalization of the bare coupling 
does not appear as a way to absorb any divergence in the self-consistent mass but as a way to avoid triviality 
(sections MI Bl and MI C[) . As it is well known, the price to pay for defining a non-trivial scalar theory by means of 
coupling renormalization is that the theory needs to be regarded as a cut-off theory, with the cut-off A taken below a 
certain scale, known as the Landau scale A p . It is then important to discuss how calculations can be made insensitive 
to A in this context. This is discussed in section MIDI where we provide a complete analytical discussion of the 
solutions of the gap equation for any value of A, below and above the Landau scale. Finally, in section MI El we 
present a new look at the gap equation based on "evolution" or "flow" equations for the thermal and field dependence 
of the self-consistent mass: these equations not only provide an efficient way to solve the gap equation but they also 
shed a new light on its renormalization. Section IIVI is then devoted to the analysis of the effective potential in the 
Hartree approximation. Its renormalization is described in sections IIV Al and IIV Bl with an emphasis on the links with 
the more general approach presented in [35| and the distinction between renormalization and consistency conditions. 
Sections IIV CI and IIVDI discuss analytically the shape of the effective potential as the temperature is lowered from an 
initial "high" temperature T+ where the theory is chosen to be in the symmetric phase, down to T = 0. We prove 
in particular that, depending on the values of the parameters, there is either no transition or a first order phase 
transition. Section [V] is devoted to a discussion concerning the renormalization in the broken symmetry phase and to 
conclusions. 
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II. THE 2PI-RESUMMED EFFECTIVE POTENTIAL 

In what follows, we consider a real scalar ip 4 theory in four dimensions at finite temperature, defined by the 
Euclidean action: 

S[<P] = I''* dr J d* X {^-(d^f + \{^) 2 + !f <p 2 + % ^ , (1) 

where the inverse temperature sets the size of the compact interval for temporal integration, and mo and Ao denote 
respectively the bare mass and the bare coupling. Notice that, in order that the spectrum be bounded from below, 
one should restrict to Ao > 0. In what follows we restrict our attention to the case Ao > 0. 

The two-particle-irreduciblc (2PI) formalism provides a representation of the effective potential "f(<fi) corresponding 
to S[ip] in terms of 2PI diagrams. More precisely, it is obtained as the value taken by the 2PI functional: 

7 [0, G] = ^ <j> 2 + ^ 4 + i [ T In G- 1 + J [ T (Q 2 + m 2 ) G + $[</>, G; A ] (2) 



(3) 



2 4! 2 J Q 

at its stationary point G = 6$, that is 7(0) = "/[(/), G^] with 
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In Eq. ([5]), </> represents a homogeneous field configuration and G(icj n ,q) a function of the Matsubara frequency 
u>„ = 2-ktiT and the three dimensional momentum q. We have also adopted the notation: 

[ T f(Q) = T Y2 I f(i"n,q)=Tj^ ! 7^ /(^n,9). (4) 

JQ „ Jq n J 

Finally, the functional — $[0, G; Ao] corresponds to all 0-leg 2PI diagrams that one can draw in the shifted theory 
S[4> + if] — S[4>] — (SS/S(j))ip at finite temperature with propagator G. This functional cannot be computed exactly. 
So-called ^-derivable approximations consist in retaining in G; Ao] only certain classes of diagrams. In this paper 
we consider the well known Hartrec approximation which corresponds to the truncation: 



T 
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$[0,G;A o ] = ^0 2 / G+^( 



According to the above discussion, in order to compute the corresponding effective potential, we first need to determine 
the stationary propagator G&. 



III. GAP EQUATION 



The stationarity condition (j3|) can be expressed equivalcntly as 

2<5$ 



G^t 



(Q) = Q* + - 



8G{Q) 



(6) 



where Q 2 = u> 2 + q 2 . We have used the subscripts <fi and T to stress the fact that the solution G^ y t depends on both 
the field <j> and the temperature T. In what follows, we shall omit this notation, unless specifically needed. In the 
case of the Hartree approximation ([3]), the propagator takes a very simple form, namely G(Q) = 1/(Q 2 + M 2 ), with 
the mass M obeying the so-called gap equation: 



AI 2 =ml + ^ 



G 



(7) 



We shall explain how to solve this equation later, in Section fill El Before doing so, a fundamental difficulty needs to 
be bypassed, namely the fact that the gap equation only makes sense within a given ultraviolet regularization. It 
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follows that its solution(s) can depend strongly on the chosen regularization, unless something is done to remove this 
sensitivity. This is what renormalization is all about in this context. 



Let us choose for instance a three dimensional, rotation invariant regularization. After performing the Matsubara 
sum, the sum-integral entering the gap equation reads explicitly: 1 



)<A Jq<A 



2e„ 



where e q = (q 2 + M 2 ) 1 ! 2 and n(e) = \j(eP e — 1). The term corresponding to "1" in the numerator of Eq. 
particularly sensitive to A. An explicit calculation leads to 



1 

<j<A 2£<j 



1 

8^ 



A V A 2 + M 2 - M 2 Arcsinh 



(8) 



is 



(9) 



which shows that the gap equation ([7]) possesses terms which depend quadratically and logarithmically on the scale 
A. What is meant by this is that if we would take A to infinity, for fixed m 2 ,, Xq, and M 2 , the right-hand-side of the 
gap equation would diverge quadratically and logarithmically. However, one should keep in mind that not all these 
divergences lead to a divergence of the solution of the gap equation. 2 The purpose of the remainder of this section is 
to clarify the connection between the sensitivity of the gap equation to the scale A and the sensitivity of its solution(s) 
to this very same scale, and to explain how the later can be removed, or at least substantially reduced. 



A. Quadratic divergence 

Let us first prove that the quadratic divergence of the gap equation ([7]) leads to a divergence of its solution M 2 . 
We first prove that there is a unique solution for large enough A and then that this solution diverges with increasing 
A. To this purpose, it is convenient to write the gap equation as = /a(M 2 ) with 



f A {M 2 ) = -AI 2 + m 2 + 



A 


2 + 


f G l 


2 




JQ<A 



(10) 



and G = 1/(Q 2 + M 2 ), and study the positive 3 zeros of f A (M 2 ). The derivative of this function with respect to M 2 
reads 



I'a(M 2 ) = -1 £ f G 2 , 

1 JQ<A 



(11) 



which is always strictly negative because Aq > 0. It follows that /a (AT 2 ) decreases strictly from 



/A(0)=m 2 + ^ 



Q<A 



(12) 



to /a(oo) = — co (the regularized tadpole integral in Eq. (jTU)) is suppressed for M 2 3> A 2 , T 2 ). Then, the existence of 
a solution of the gap equation depends on the sign of /a(0). Even though the parameter m 2 , could be negative, the 
quadratic and positive divergence in /a(0), given explicitly in Eq. ©, implies that there exists a value of A above 
which /a(0) > and the gap equation admits a unique solution M 2 . For the same reason, given a mass [i, there 
exists a value of A above which f\ (fj, 2 ) > and thus M 2 > fi 2 . This shows that the solution M 2 diverges as A — > oo, 
for fixed m 2 , and Aq, as announced. 



1 Depending on the context, we use the notations Q < A and q < A to designate the same three dimensional regularization. 

2 It is true that these divergent terms lead to quadratic and logarithmic divergences in the coefficients of the formal pcrturbativc expansion 
of the solution M 2 in powers of Aq. However it is not true that they all lead to a divergence of M 2 , see the subsequent discussion. 

3 We shall only be concerned with positive (M 2 > 0) solutions of the gap equation. For a discussion of negative solutions of the "explicit" 
form of the gap equation at zero temperature, that is the gap equation where the integrals have been performed explicitly for M 2 > 
and then extended to any value of M 2 , see [52H . 
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B. Renormalization of the mass — Triviality 

Since the quadratic divergence of the gap equation depends neither on <j) nor on T, it can be absorbed by adjusting 
the divergent part of TOq. On the other hand, the finite part of m 2 , can be used to impose a condition at a given value 
of <p an d a given value of T. We choose (f> = 0, T = and impose the renormalization condition: 



Mi 



=0, T„ 



(13) 



The parameter m 2 is positive by construction, because it is a solution of the gap equation at <j> = and T = T+. We 
shall choose it strictly positive in what follows. The renormalization condition (|13[) is quite natural when studying 
how the system evolves as the temperature T is decreased from a "high" temperature T+ where the system is required 
to be in the symmetric phase, see Section llVl It can be rewritten as a choice of the bare mass, namely: 



Ao 
2 



T, 



G* 



Q<A 



(14) 



with G* = 1/(Q 2 + m 2 ). With this choice, the gap equation reads = f\(M 2 ), with 



/a(M 2 ) = -A/ 2 + m 2 + 



Ao 
2 



G — G* 

)<A JQ<A 



(15) 



The dependence of the gap equation on A has been changed by the renormalization procedure and we need to 
reconsider its possible solutions and their dependence on A. As before, we do so by discussing the zeros of the 
function /a (A/ 2 ). 

Notice first that f' A (M 2 ) = f' A (M 2 ) is strictly negative in view of Eq. ([IT]). It follows that the function f\(M 2 ) 
decreases strictly from 



/a(0) 



Ao 
2 



Q<A 



Q 2 



<A 



(16) 



to /a(oo) = — oo (the regularized tadpole integral in Eq. (fT5")) is suppressed for M 2 A 2 , T 2 ). Using Eqs. (JS) and 
©, it is easily checked that /a(0) diverges logarithmically like cm 2 In A, with c > 0. Then, there exists a value 
of A above which /a(0) > and the gap equation admits a unique solution Af 2 . For the same reason, for any 
Am 2 > 0, /a(^ 2 ± Am 2 ) diverges logarithmically like =FcAm 2 lnA. Then, there exists a value of A above which 



/ A (m 2 - Am 2 ) > and f\(mi + Am 2 ) < 0, and thus \M 1 



for fixed Ao- Using the gap equation 
this limit is approached. We obtain 



/a(M 2 ) as well as Eqs. 



< Am 2 . This proves that M 2 — > m 2 as A — > oo, 
|) and (JUJ) , it is possible to determine precisely how 



M 



8tt 2 


2 + 




In A 




J q £ q . 



(17) 



where e* = (q 2 + m 2 ) 1 / 2 and 8*n e = n e — n*, with n* the thermal factor at temperature T±. Thus, unlike what 
happened with the quadratic divergence, the logarithmic divergence of the gap equation (the one that remains after 
mass renormalization) does not lead to a divergence of the solution M 2 . Rather, if we insist in taking the limit 
A — > oo for fixed Ao, the solution of the gap equation converges to the renormalized mass m 2 for any value of the field 
(f> and the temperature T. This illustrates the triviality of (p 4 theory (5lj . at least in the particular approximation 
considered here. Equation (|1T[) shows that the trivial limit is approached rather slowly (logarithmically). 



4 The absence of divergence in the solution of the mass renormalized gap equation means that the resummation of supcrdaisy diagrams 
encoded in this equation formally resums the perturbative logarithmic divergences associated with each diagram of this series into a 
convergent contribution. That perturbative divergences can be resummed to non-divergent expressions was observed in I42l l53ll . 



C. Renormalization of the coupling — Landau pole 



The previous analysis has shown that, for the particular choice of mass renormalization we have considered, the 
triviality of <p A theory is related to the presence of a logarithmically divergent term in the gap equation. In order to 
define a non-trivial theory, one can absorb this divergence in a redefinition of A . Using the formula (sec Appendix [A"l): 



as well as the identity: 



G = 



Q<A 



G 



<A 



q<A £ q 



(18) 



G - G + = -{M 2 - ml) G* G = -(M 2 - ml) G 2 + (M 2 - m 2 ) 2 G 2 G , 



(19) 



it is clear that the logarithmic divergence in Eq. (fT5j) is entirely accounted for by the term —(A/ 2 — m 2 ) Jq* <a G 2 
Isolating this contribution and defining an effective coupling A* [54| such that 



1 



i 

To 



G 2 , 



Q<A 



(20) 



the gap equation = Ja(M 2 ) can be rewritten as = g\(M 2 ) with 



g A (M 2 ) EE -M 2 + ml + ± 



-M 



A* 
2 



G - 

<A 
q<A £ q 



G + + (M 2 - ml) 



Q<A 



+ (M 2 -m 2 ) 2 



T, 



G 



<A 



G , G 



Q<A 



(21) 



where, in the second line, we have made use of Eqs. (118)) and (fT9|) to obtain an explicitly convergent expression that 
we will use later. If we insist in keeping the bare coupling Ao fixed and positive, then, according to Eq. (|20[) . the 
effective coupling A* goes to zero as A — > oo and we recover the trivial result with M 2 — m 2 ~ A* [tp 2 + f <5*n e * /e*]/2, 
in agreement with Eq. (JTTj) . In contrast, fixing the value of A* allows to avoid triviality. This corresponds to the 
following redefinition of Ao (Ao is then <p- and T-independent as it should): 



1 

Ao 



A, 



G 



(22) 



Q<A 



However, it appears that in order to maintain Ao > 0, and in turn a meaningful microscopic theory, one needs to 
restrict the cut-off A below a certain scale known as the Landau scale or Landau pole A p , defined by 



G 



(23) 



Q<A P 



In other words, the non-trivial Hartree approximation has a meaning only if it is considered as describing an effective 
theory. Notice that, from Eq. (f2"3"]). it follows that A* > (negative values of the renormalized coupling could 
be possible with other renormalization schemes without violating the requirement Ao > 0). More explicitly, using 
Eq. (fMJ), we obtain 



1 

aT 



i 



16tH 



Arcsinh 



A, 



A, 



e q n 



q<A p 



2e* 



(24) 



In the limit A* — > + , A p — > oo with A 2 ~ /i 2 e 327r / A * where (i* depends on m* and T*: 
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D. Ultraviolet sensitivity 

Because, in the Hartree approximation at least, a non-trivial ip 4 theory is only to be considered as a cut-off theory, 
we need to wonder how calculations done within such a theory can be made almost independent of the cut-off A. 
The main idea is that, if one considers the regime A p 3> T+, m*, <fi, T, one can choose the cut-off A such that both 
requirements A p > A and A 3> T*, m*, <f>, T are met. Then, because all divergences of the gap equation have been 
absorbed in the redefinition of the bare parameters, we expect that its solution M 2 be almost cut-off independent, 
up to terms of order T*/A , to*/A , 0/A , T/A <C 1. The Hartree approximation offers the possibility to illustrate 
this issue, for one can study its solutions as a function of A, in particular as A — > oo, and thus assess under which 
conditions these solutions can be considered practically insensitive to A for A < A p . 



4>p(T,' 



4>c("0 












1 






2 







1 



A=oo 



FIG. 

A, Ac 



1: Number of solutions of the gap equation = g\{M 2 ) in the (A, 2 )-plane in the regime of interest, that is when 
S> T*, m*, 4>, T. Orange areas correspond to no solution, green areas correspond to one solution, and the blue area 



corresponds to two solutions, as indicated also by the labels. In the regime of interest <j(- 2 ~ A 2 /(8tt 2 ) and <j>p(T, oo) ~ A 2 /(47r 2 e 2 ) 
are both high scales. 



The analysis of the solutions of the gap equation = g\(M 2 ) as one varies the cut-off A is given in Appendix [B] 
where, in order to simplify the discussion, we restrict our analysis to the "regime of interest" that is a regime where 
both A and A p are much larger than all the other scales T*, to*, <f>, and T. The result of this analysis is that the number 
of solutions of the gap equation depends on the values of A and cj> 2 , as we represent in Fig.[T]and explain in what follows: 

1. For A < A p , there is a "critical" value cf> 2 (T,A) of <fi 2 such that the gap equation admits a unique solution if 
4> 2 > 't'liT; A) and no solution if <f> 2 < (j) 2 (T, A), see Fig.[T] Notice that <fi 2 (T, A) is not necessarily positive. In practice 
it is useful to know when it could become strictly positive, signaling the fact that the gap equation has no (positive) 
solutions if < (f> 2 < <fi 2 (T,A). To this purpose, one writes, see Appendix |B| 

£(r,A) = -^- / 5a, (26) 

A * Jq<A 1 

with 

C* - ml + ± 



q<A 9 



Q<A 



Q 2 



(27) 
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Clearly, (f> 2 (T,A) is strictly positive if the parameters A, T+, m*, and A* are such that C\ < and the temperature T 
is strictly below a certain "critical" temperature T c defined by 



2 a 



(28) 



Using Eq. (f2"T|). it is easy to check that </> 2 (T*, A) is strictly negative. It follows that T c < T+. Moreover, in the regime 
of interest, one can neglect the dependence of C* and of the thermal integrals with respect to A. The critical value 
<f% (T, A) is then independent of A (this is the reason why it is represented by an horizontal line in Fig. [T]) and one 
has </> 2 (T) = (T 2 - T 2 )/12 with T 2 = -24C*/A*. 

2. For A = A p , a new "relevant" scale appears: 



G* >0, 



(29) 



<A D 



see Appendix IB"1 It is such that (f> 2 > </> 2 (T) an d the gap equation admits no solution if <fi 2 > (f> 2 , one solution 

< (f> 2 (T), see Fig. [T] Notice that, in the regime of interest, 



if 



> 2 > 0c (-^)> an d again no solution if 
6 2 ~ Ap/(87r 2 ) is a high scale. 



3. Similarly, for A > A p a new relevant scale (pp(T, A) > <fi 2 (T) enters the discussion. The gap equation admits 
no solution if 4> 2 > 4> 2 (T,A) one solution if 4> 2 = </> 2 (T, A), two solutions if 2 (T, A) > 2 > </> 2 (T) and one solution 
if 4> 2 < 4> 2 (T), see Fig. [T] The value of (j) 2 (T, A) as A -> A+ is nothing but 2 . Moreover, in the regime of inter- 
est, one can show that <pp(T, A) decreases with A, from <fip(T, A+) = <fi 2 to <pp(T, oo) ~ A 2 / (47r 2 e 2 ), another high scale. 

The previous analysis shows in particular that the unique solution of the gap equation which exists for A < A p and 



if) 2 > 4> 2 (T) and which we shall call the physical solution, behaves rather differently, depending on the value of <f. 
one increases the value of A. If <fi 2 > (f> 2 as we show in Appendix [Bj the physical solution diverges as A — > A p 



as 
In 



contrast, if cf> 2 > <f> 2 > <f) 2 (T), nothing of this kind occurs: the physical solution can be extended beyond the Landau 
scale and a second solution, which we call the unphysical solution, appears. If 2 > (f> 2 > </> 2 (oo), there is value of A at 
which the two solutions merge into one solution and then cease to exist (without diverging). On the contrary, in the 
convergence window (f> 2 (oo) > tfi 2 > <fi 2 (T) the two solutions can be extended to A — > oo and one can then assess how 



rapidly they converge to their limiting values 
with respect to A. We obtain 

DM 2 



To do so, we take a derivative of the gap equation = g\{M ) 



dgA 



dA 



g' A (M 2 ) dA 



(30) 



where dg\/dA corresponds to the explicit dependence of g\ with respect to A. An implicit dependence is also present, 
through M 2 , which explains the appearance of g^(M 2 ) in the previous formula. We are allowed to divide by g' A (M 2 ) 
because, as long as <fi 2 < <fi 2 (oo), one has g' A (M 2 ) ^ 0, see Appendix |B] Moreover, this remain true in the limit 
A — > oo. It follows that 



dM 2 
~~dA 



1 



_ dgA 



A/2 = Ml 



Then, a simple calculation using Eqs. (|AT|) . (|A"4]) and ||BlJ leads to 



dM 2 
~~dA 



647rVoo(^ 



A :i 



(31) 



(32) 



or, in other words: 



M 2 
m 2 



M 2 



3A* 



1287r 2 |^(AO 



A 2 ' 



(33) 



where we have expressed all the masses in units of the renormalizcd mass to*. Although we did not make it explicit, 
the limiting value depends on T and <f>. In the case of the physical solution increases strictly from at 
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4> 2 = 4> 2 C {T) to a very large value Mp(oo) ~ 4A 2 /e 2 at <f> 2 = 2 (oo) - A 2 /(47r 2 e 2 ), see Appendix |B] Thus, in the 
regime of interest, such as in particular <f> <C A p , M^/m 2 is of the order 1 or smaller and the scale multiplying the 
convergence factor 1/A 2 in Eq. (p3|) is a low scale of the order of m 2 <C A 2 or smaller (the factor l/g'^M^) can only 
improve the convegence, for it goes to zero when — > 0). Then, the insensitivity of the physical solution to the 
scale A is expected to be observed already for values of A below the Landau scale A p . This expectation is confirmed 
by the numerical results presented in Fig. [2] which shows the convergence of the physical solution as A — > oo for 
increasing values of <fi '■ for small values of the field the physical solution has already converged for values of A well 
below the Landau pole, while for very large values of <f> the convergence, if it happens, occurs for values of A above 
the Landau pole. This is a clear illustration of why, although there exists a Landau pole A p which forces us to 
choose the cut-off A < A p , one can define, in some region of the parameter space, a theory whose results are almost 
insensitive to the cut-off A. 




FIG. 2: Convergence of the physical solution of the renormalized and A-dependent gap equation at different values of the field 
for ml/T? = 0.5, A* = 20 and T/T* = 0.012. The dahed vertical line, common to both plots, indicates the location of the 
Landau scale. The left panel shows the dependence of the rescaled physical solution with respect to the cut-off A. The right 
panel illustrates the validity of the asymptotic formula (|33p represented by the solid lines. The inset shows the behavior of the 
solution for two values of the field which lie outside of the convergence window (see text for details). 



E. Flow approach to the gap equation 



To conclude this section, we present an elegant way to solve the gap equation. This approach is also interesting 
on formal grounds because it gives a new perspective on renormalization. It could also be used for higher order 
approximations. 

Let us assume that the solution of the gap equation is known for some 4» an d T, and let us ask how this solution 
changes as we vary <f> or T. The change in the solution is governed by "evolution" or "flow" equations that we now 
derive. Let us first express the infinitesimal change of M 2 with respect to tf> 2 . From the original bare gap equation 
([7]) and using the fact that the bare parameters m 2 , and Ap arc (^-independent, we obtain 



dM 2 



Ao 
2 



T 
Q<A 



dG 



= g 2 dAP 



Q<A 



(34) 



where we used SG 
rewritten as 



-G 2 5G 1 . Collecting the derivatives of M 2 on the same side of the equation, this can be 



dM 2 



V 
2"' 



(35) 
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with the function V defined by 

1 1 1 r T 

G 2 . (36) 



V A 2 



<A 



Similar manipulations are possible to express the infinitesimal change of M 2 with respect to T. Using the fact that 
the bare parameters are T- independent, one writes 

dM 2 A / d f T \ A f T ^ dM 2 

~~dT~ ~ ~2~ ydT Jq < a J 2~ Jq < \ ~W> (37) 

where the notation (d/dT Jq <a )G defined in Appendix |A1 refers to a derivative with respect to the explicit thermal 
dependence, not the implicit one, hidden in G. As before, this can be rewritten as: 



dM 2 V ( d f T 

Iq<a. 



dT 2 \ dT 



G . (38) 



A similar strategy can be used to express the infinitesimal variations of V itself. We obtain: 



G 3 , (39) 




dV 

The interesting feature of Eqs. ([35]) . (|38|) . (|39|) and (|40|) is that they are explicitly finite, without the need to impose any 
renormalization condition. They can be used as an alternative tool to solve the gap equation. The only information 
that remains to be specified are the initial conditions for M 2 and V, for instance at <f> 2 = and T = T*. The initial 
condition for M 2 reads M?_ T< = m 2 which is nothing but the renormalization condition (|13p . As for V, its initial 
condition reads: 

fy=o,T 4 = A* , (41) 

and it is easily checked that it can be reinterpreted as the the renormalization condition corresponding to 
the choice of Ao in Eq. (|22|) . Thus, specifying the initial conditions in the flow approach corresponds to impos- 
ing renormalization conditions in the standard approach 5 . As already discussed, m 2 and A* are taken strictly positive. 



In practice, it is more convenient to solve the evolution equations (|35l) and (|38[) coupled to the explicit and renor- 
malized expression for V which we obtain from combining Eqs. (|36|) and (|22[) : 



1 1 1 

v ~ a7 + 2 



T f T„ 

G 2 - Gl 

Q<A J Q<A 



(42) 



Because 4> 2 (T+) < 0, we first solve the equation (|35|) at T = from the initial condition M| =0 T t = m * > ^ and 
obtain the function Mj Tt . For each value of <j>, the solution M| T ^ is then used as the initial condition for the flow 
equation (|3"8")l . which we solve from T* for decreasing values of T. We obtain in this way the function Mj T for any 
value of (j) and T (with </> 2 > 4> 2 (T)). This method will be used in the next section where we evaluate the effective 
potential. Notice finally that V = — X* / g' A (M 2 ) where g\(M 2 ) is the function introduced in Eq. (|2ip to discuss the 
solutions of the rcnormalized gap equation (see Eq. (|B2I) for its derivative). We are only interested in the physical 
solution to the gap equation. Irrespectively of the value of A, this solution is such that g' A (M 2 ) < 0, see Appendix IB1 
and thus V > 0. From Eqs. ([33)) and (|3"9")l , it follows then that M 2 and V increase with (j) 2 . In particular at T = T + 
and for any (j) 2 , one has M 2 > m 2 and V > A*. We will make use of this remark later. 



5 This is similar in spirit to the approach followed in [HEIl . 
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IV. EFFECTIVE POTENTIAL 



Once the mass M and the propagator G have been determined, one can evaluate numerically the renormalized 
version of the effective potential to be given in Section IIVBI and study its change of shape as the temperature T is 
lowered from the "high" temperature possibly triggering a change of phase. In this section we show that a large 
part of the analysis can be done analytically. We prove in particular that in the Hartree approximation, for those 
values of the parameters such that there is a phase transition, the later cannot be of second order. To this purpose, we 
use information on the geometry of the effective potential, such as its curvature, as well as information on the location 
of its possible extrema, as given by the field equation. Notice that our analysis will be performed in the presence of a 
cut-off A below the Landau scale A p . After renormalization of the effective potential, see below, and for parameters 
such that A p is much larger than all other scales in the problem, the renormalized effective potential computed with 
A < A p will be essentially the same as the one computed in the limit A — > oo. We shall use this simplification when 
evaluating the effective potential numerically. 



A. Geometry of the effective potential 



The geometry of the effective potential is encoded in its field derivatives. For instance, the curvature for a given 
value of 6 reads: 



M 2 S = ml 



X 



0,2, ^0 



G 



An 



<A 



T G 2 ^ 

KA 



d<j) 



(43) 



where we made use of the stationnarity condition ([3]) and 5G = —G 2 5G 1 . From the flow equation (|35[) . one obtains 
DM 2 /d<j) = V<f>. Plugging this into Eq. (|4"3")l and using Eq. (|3l)|). we arrive at: 



M 2 = M 2 + V <f) 2 - An 2 . 
Similarly, the fourth derivative of the effective potential at = is given by: 



v ^ 



The same argument as for M 2 leads to: 



= Ao - - Ao / G^ =0 

6=0 Z JQ<A 



d 2 M 2 



4=0 



Vo 



<f>=0 



SV a 



= 



2A . 



(44) 



(45) 



(46) 



These results illustrate the general discussion given in the Introduction concerning the existence of different definitions 
of n-point functions in a given approximation of the 2PI effective action. For instance, the two-point function (here 
simply a mass) could be defined cither from the solution M of the gap equation or from the second derivative M 2 
of the effective potential. Although these two definitions should coincide for any value of an d T in the absence of 
approximations, Eq. (|44|) clearly shows that, in the Hartree approximation, they coincide for = only. For other 
approximations, the discrepancy between the different definitions of the two-point function can also be observed for 
= 0. Similar remarks apply to the four-point function. Notice however that these discrepancies always appear as 
higher order effects. In the present case for instance, if we consider a formal expansion in power of the bare coupling 
Ao, we obtain M 2 — M 2 = O(Xq) as well as V — V = O(Xq) which show that the discrepancies are beyond the 
accuracy of the Hartree approximation. 



More generally the discrepancies between the different definitions of a given n-point function are always of higher 
order as compared to the accuracy of the 2PI truncation that one considers and are thus of little relevance a priori. 
However, this is only true for n > 6. For n = 2 and n = 4, the discrepancies between the different definitions of 
n-point functions can be divergent, and they are thus negligeable only after divergences have been properly absorbed. 
For instance, the two different definitions of the two-point function can bring independent quadratic divergences. In 
order to remove them, one needs to allow for two independent bare masses following the general prescription of [35| . 
It turns out that in the case of the Hartree approximation, the quadratic divergences are not independent but equal 6 



This has to do with the fact that the two definitions of the mass coincide for rf> = 0. 
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which allows one to work with only one bare mass. Nevertheless the sensitivities M 2 and M 2 with respect to A are 
radically different for, as we have seen above, M 2 is defined for any value of A with A-dependent terms of the form 
1/A 2 , whereas, due to the presence of Ao in Eq. ([44]) . M 2 diverges as A — > A~. For the same reason, the sensitivities 
of V and V with respect to A are different. One possible way out is to argue, in a way similar to (46|, that if one 
considers A< A p , but much larger than all other scales in the problem, the additional sensitivity of M 2 and V as 
compared to M 2 and V is a very mild one. Here we follow a different approach, 7 in line with the general presentation 
in [35j . which is generalizable to higher order approximations. Since the difference of sensitivity with respect to A 
comes from the fact that the two definitions V and V of the four-point function differ within a given approximation, 
it is natural to introduce a bare coupling A4 different from Aq in such a way that the original 2PI functional reads 



\ f InG^ + J f (Q 2 +m 2 )G 
1 Jq<a 1 Jq<a 



^0 ,9 f T , A 



4 




f G + ^| / G ) . (47) 

<A 



Clearly Ao enters the definition of V, as before, whereas A4 appears as a tree level contribution to V. With this 
modification, using the same steps as those leading to Eqs. (j4"4"|) and ([41)]) we obtain 

M 2 = M 2 + V4> 2 + i(A 4 -3A o )0 2 , (48) 

V^=o = 3F 0= o + (A 4 - 3A ) . (49) 

We can now adjust A4 such that V is as insensitive to A as V, that is up to terms of order 1/A 2 . This however does 
not fix the finite part of A4. In order to fix the later without introducing additional parameters, and following [35ll44|. 
we impose the consistency condition: 

fy=0,T. =^=0,T., (50) 

at the renormahzation point = and T = T*, condition which should hold automatically if no approximation was 
considered. This condition translates into a particular choice of the bare parameter A4: 

A 4 - 3A = -2A* . (51) 

Using this particular choice, the derivatives of the effective potential become then 

M 2 = M 2 + (V — X*)4> 2 , (52) 
V^ Q = 3t^= -2A*. (53) 

Notice that M 2 and V are now both cut-off insensitive: the renormahzation of V automatically took into account 
the renormahzation of M 2 . We stress that this is only specific to certain approximations. In general, one would need 
to introduce a second bare mass which one would fix by imposing another consistency condition. Notice also that, 
if we expand the bare couplings in powers of the renormalized coupling A* using Eqs. (|22p and (|51[) , we find that 
A4 — Ao = C(A 2 ), in line with the fact that the artifacts that the introduction of a second bare coupling A4 is meant 
to cure lie beyond the accuracy of the Hartree approximation. As a final remark, recall that, as we have shown in the 
previous section, M 2 > ml > and V > A* > at T = and for any 4> 2 . It follows that M 2 > ml > at T = 
and for any <fi 2 . This means that the effective potential at the initial temperature T = T* is strictly convex. 

B. Renormalized effective potential 

We are now in a position to derive a finite expression for the effective potential. Using Eq. (JT)) in Eq. (|47|) the 
effective potential 7(0) = j[cf>, G] reads 

7 (0) = i(A 4 -3Ao)0 4 + ^T^ + ^r k ^4r Q2 °- ^ 
41 ^A I J Q<A I J Q<A 



7 Despite this different point of view, our conclusions regarding the order of the phase transition will be the same as those in [4&I . 
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Using Eqs. (|36|) and (|5Tj) . we arrive at 



7(0) 



2A^ ~ 12 



M 4 
2U 



Q<A 



In G- 1 + Q 2 G- -M 4 G 2 



(55) 



Expanding the integrand for large Q we observe that it has no quadratic or logarithmic divergence, only a quartic 
one of the form /<p <A (l n Q 2 + !)• We can thus write 



A* 



M 4 



rw 2A 12 v 2F 2 



<A 



ln(Q J + M ) - In 



M 2 G M 4 G 2 

2 



7' 
Q<A 



lnQ 2 



1 



(56) 



The first integral is convergent. The second is divergent, but one can convince oneself that its divergence is T- 
and (^-independent. Then, if we subtract the effective potential evaluated at T+ and (f> = 0, we obtain an explicitly 
convergent expression 



7(^-7^(0) = -^ 



A* ,4 



M 4 

w 



2K 



q<A 



q<A 



q<A 



£ q -q- 



2Ta 



-2Tln(l - e 
-2T*ln(l 



M 4 ' 


1 




~ 2 








+ M 2 



e - £ 9 / T *) + m \ 



q<A 

fc q 
Sq 



£ 



M 



8ef 



4^3 



e* n*' , 



4e* 



(57) 



which is our final form for the renormalized effective potential. 8 This formula for the effective potential differs from 
that in (4(| in that we have used the additional bare parameter A4. The approach we have followed here to obtain a 
finite effective potential can be generalized to higher orders of approximation, see for instance [40]. 



C. Extrema of the effective potential 

At any temperature, the extrema of the effective potential (I47|) are given by the field equation: 

0=P 



(58) 



In making this equation explicit one can take advantage of the stationarity condition ([3]) . In the Hartree approximation, 
and taking into account the newly introduced bare coupling A4, we arrive at: 

°^( m « + lP + Y /J <A ^) ' (59) 

Using the relation (|5"Tj) as well as the gap equation (0, we finally arrive at the system of equations: 



Q = 4>[-fP -MfJ and = ffA (M|,^), (60) 

where we have made explicit the dependence of g^(M 2 ) with respect to <j>. Notice that this system of equations is 
explicitly finite. This is not surprising for, once the second and fourth derivatives of the effective potential have been 
renormalized, the effective potential is renormalized (up to an overall constant), and so is any information that can 
be extracted from it such as the location of its extrema. 9 



8 The first two integrals can be done analytically with the result (M 4 — m 4 ) / (128n 2 ) in the limit A — > 00. 

9 Notice that the way we have obtained an explicitly renormalized field equation differs in spirit from that followed in l5rj| for the case 
of the O(N) model. There, the field equation was similar to the one in Eq. 16011 with A* replaced by the bare coupling (denoted by 
A in Ref. [53] )• Since the latter had no divergence to absorb in the field equation, it was then identified to the renormalized coupling 
(although this identification was not possible at the level of the gap equation). Even though this result is correct and agrees with 
Eq. l |60|l . its justification in [5t3 | hides the origin of the difficulty, namely the existence of multiple definitions for n-point functions and 
does not allow a generalization to higher orders of approximation. For a discussion of similar issues, see also [57| . 
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According to Eq. (|60)) . the trivial extremum = only makes sense as long as the gap equation at (f> = has a 
solution. This ceases to be true if </> 2 (T) is strictly positive, that is when C± < and T < T c , see our earlier discussion. 
On the other hand, the non trivial extrema obey the system: 



which we rewrite conveniently as 



0=y^ 2 -Mj and = 5A (M|> 2 ), (61) 



and = /i A (M|), (62) 



where the function h A (M 2 ) does not depend on the field and is given by 

_3_ 



h A {M' 2 ) = g A [ M 2 , ^-M 2 ) = |m 2 + g A {M 2 , 2 = 0) , (63) 



with first and second derivatives equal to 

h' A (M 2 ) = l+g' A (M 2 ), (64) 

h'L{M 2 ) = g'i(M 2 ). (65) 

The non trivial extrema of the effective potential can thus be obtained from the zeros of h A (M 2 ), which we discuss 
below. Notice that, from Eq. (|52|) . the curvature of the effective potential at the non-trivial extrema is given by 



2K\ - 2 . / 1 2\ - 2 2Kh' A (M 2 ) 



Because g' A (M 2 ) < when A < A p , see Appendix [Bl the sign of the curvature is the sign of h' A {M 2 ). 



D. Temperature dependence of the effective potential 

Let us now discuss the zeros of h A (M 2 ) and the corresponding curvature encoded in Eq. (|66j) and deduce the shape 
of the effective potential as we decrease the temperature from T* down to T = 0. Because h' A (M 2 ) > 0, the function 
h' A (M 2 ) increases strictly from h' A (0) = — oo to 

h' A (oo) = | +g' A (oc) = = ! + £ p Gl>0, (67) 

Z Z At) Z Z JQ<\ 

where we have used (?^(oo) = — A*/Ao, see Appendix IBl and Eq. (f2"2"j) . The function h' A (M 2 ) has to vanish for some 
M 2 . It follows that the function h A (M 2 ) has a minimum at M 2 = M 2 : h A (M 2 ) decreases strictly from h A (0) 
to h A (M 2 ) and then increases towards h A (oo) = oo. The number of non-trivial extrema of the effective potential 
depends thus on the signs of h A (0) and h A (M 2 ) which we now discuss. 

Let us consider h A (0) first. It is nothing but g A (0) at <j) = 0, sec Eq. (|63|) . and thus from the discussion above 
Eq. JB6|): 

M0) = ~ <fil (T, A) = a + ^ / ^ , (68) 

Z Z Jq<A Q 

where C* was defined in Eq. (|27|) and depends on the parameters. As we have seen below Eq. (|28|) . <j> 2 (T, A) is strictly 
negative for T = T+. It follows that h A (0) is strictly positive for T = T+. Moreover, the previous formula shows that 
h A (Q) decreases strictly as we lower T and reaches C* at T = 0. If the parameters are such that C* > 0, then h A (Q) 
remains strictly positive down to T = 0. If C* = 0, h A (0) remains strictly positive down to T > and vanishes at 
T = 0. If C* < 0, h A (0) remains strictly positive while T > T c (0 < T c < T*), vanishes at T = T c , changes sign and 
remains strictly negative down to T = 0. The same type of analysis needs to be done for h A (M 2 ). Using Eq. (|21l) . 
we obtain 

hA (M 2 )= 1 -M 2 + m l + ± 



8±n e 



+ (AI 2 ~ml) 2 



q<A 



(69) 



Q<A 
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FIG. 3: The temperature evolution of the effective potential obtained for ml/T 2 =0.1, A* = 3 in the limit A — > oo by evaluating 
Eq. (J57)) with M 2 calculated at T* using the ^-flow Eq. (@5) and at T / using the T-flow Eq. $38$. On each plot, a constant 
has been subtracted for convenience, without affecting the physical interpretation. 



Clearly, h\(M c ) > at T = T*. Moreover, because h' A (M e ) = 0, the thermal dependence of h(M e ) is the explicit 
one, encoded in the tadpole integral. It follows that h\{M%) decreases as one decreases T. We shall denote by its 
value at T = 0. If the parameters are such that D+ > 0, /ia(M 2 ) remains strictly positive all the way down to T = 0. 
If D+ = 0, it remains strictly positive down to T > and vanishes at T — 0. If £>* < 0, h\(M^) remains strictly 
positive while T > T s (0 < T s < T*), vanishes at T = T s , changes sign and remains strictly negative down to T = 0. 

A priori, discussing the signs of /ia(0) and h\(M^) involves nine cases, depending on the values of C* and and 
thus nine regions in parameter space. However, because C* and D* are the values of h\(0) and h^(M^) at T = 
and /ia(0) > /ia(M^ ), we have necessarily > D*, which reduces the number of regions to five: 

1. Region corresponding to C* < 0: For such parameters, h^(A4^) changes sign at T s and /ia(0) changes sign at 
T c < T s . The evolution of the effective potential as we decrease T from down to T = goes as follows: 

— For T s < T < T*, both h\(0) and h\{M^) are strictly positive and the effective potential has only a trivial 
extremum at <fi = 0. It is a minimum because M 2 = M 2 > at = 0. The effective potential is convex, see 
Fig. [Sta). 

— At T — T s , h\(0) is still strictly positive but h(M^) vanishes: in addition to the trivial minimum at <f> = 0, a 
non trivial extremum appears. According to Eq. ()66[) . its curvature is zero. The only possibility is that it is an 
inflection point of the effective potential, see Fig. Efb). 

— For T c <T < T s , /ia(0) is still strictly positive and /i(M 2 ) is strictly negative: in addition to the trivial minimum 
at <j) — 0, two non-trivial extrema are present. According to Eq. (|66[) . the non-trivial extremum closest to <p = 
is a maximum and the other one is a minimum, see Fig. [3]^). As it will become clear in the next point, it 
turns out that in this temperature interval there is a temperature T t (T c < T t < T s ) at which the minima are 
degenerate and below which the non trivial minimum becomes the absolute minimum, see Fig. EJd). 

— At T = T c , h\(0) vanishes and /ia(M 2 ) is strictly negative: the non-trivial extremum which corresponded to 
a maximum is at <j> — and merges with the trivial minimum. The curvature at <fi = is zero. Because there 
is still a non-trivial minimum at 4> ^ 0, the only possibility is that at T = T c the effective potential presents a 
maximum at (f> = 0, see Fig. G3e), with the consequence that there was a temperature T t (T c < T t < T s ) below 
which the non-trivial minimum became the absolute minimum. 

— For T < T c , the gap equation has no solution at <fi = 0: the trivial extremum disappears. The effective potential 
is defined only for <fi 2 > <^ 2 (T) and has a single non-trivial minimum, see Fig. (3Jf). 
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2. Region corresponding to C* = 0: It is just a line in parameter space. The same analysis as above applies except 
from the fact that T c = 0. The final shape of the potential is that of Fig.[3Je). 

3. Region corresponding to C* > and D+ < 0: The same analysis applies except from the fact that there is no 
T c . The final shape of the potential at T = is either that of Fig. [3jd) or that of Fig. G3c), but without evaluating 
the potential we cannot decide whether or not the non-trivial minimum became the absolute one at some temperature. 

4. Region corresponding to > and D+ = 0: It is just a line in parameter space. The same analysis applies 
except from the fact that there is no T c and T s = 0. The final shape of the potential is that of Fig. [3{b). 

5. Region corresponding to C* > and D+ > 0: There are no T c or T s . The final shape of the potential is convex, 
see Fig.EHa). 

The previous discussion shows that, if there is a phase transition, it is first order due to the finite jump in the loca- 
tion of the absolute minimum of the effective potential (see Fig. [3]). The transition occurs at T = Tt while T c and T s 
correspond to the upper and lower spinodal temperatures at which a maximum and a minimum of the potential merge. 




mfn? 



FIG. 4: Parameter space and order of the phase transition. The labels on the curves indicate the value of ln(A p /m*). We focus 
on the region ln(A p /m*) > 5 which is such our results can be considered cut-off insensitive for a cut-off A below the Landau 
scale. Note that this region lies below the dot-dashed curve that is in the region of validity of the condition (|B12I) which played 
a role in part of our analysis of the gap equation. 



These conclusions are valid for a given value of the cut-off A < A p . As we have explained in the previous section, if 
we want the results to be independent of the cut-off A we need to choose the parameters such that A p 3> jto*, T*, T 
so that it is possible to choose A below the Landau scale but larger with respect to all other scales. 10 In what follows 
we choose T* = 1 which sets the mass unit. Since we are only interested in values of T below T+, it is then enough to 
choose and A* such that Ap/m^ 3> 1. For definiteness, we choose to work in the region of parameter space such 
that ln(A p /m+) > 5. Note that, in this regime 



For practical purposes, it is then enough to compute the effective potential for A = oo without introducing any substantial difference 
as compared to the effective potential at A < A p . One can also perform an analytical study of the effective potential for A > A p . The 
only subtlety is that g'(M 2 ) can be positive for M 2 > M 2 and the discussion around Eq. J66H needs to be revisited. One can convince 
oneself that the zeros of h(M 2 ) are all low scales, lying below the high scale M 2 , and correspond thus to negative values of g'(M 2 ). 
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which wc obtained from the equation h' A (M 2 ) = at T = and for A ^> m*, T+. One can show in the same way that 
which is defined as h\(M 2 ) at T = 0, is given by 

Since regions 2., 3. and 4. are such that C* > and < 0, it follows that < C* < A*M|(T = 0)/(32tt 2 ) and they 
form a very narrow band, indistinguishable from the line C* = whose equation is 11 

-l 

(72) 

The parameter space is then essentially divided in two regions, see Fig. HI The lower (grey) region, corresponding 
to region 5. of the analysis above, is such that no phase transition occurs between T = T* and T = 0. The 
upper (white) region, corresponding to region 1. of the analysis above, is such that a first order phase transition 
occurs at T = T t . The temperatures T t and T s can only be accessed numerically. We find that for moderate 
values of the parameters A* and m 2 the first order phase transition is rather weak, that is for reasonable values 
of the parameters the difference between T c and T s is small and also for T c < T < T s the non-trivial minimum 
is close to the trivial one. The first order nature of the phase transition is an artifact of the Hartree approxi- 
mation which corresponds to a 0(A) truncation of the 2PI functional at skeleton level. Since this effect is rather 
weak, an improvement of the truncation by including the order-A 2 and field-dependent skeleton contribution to 
the 2PI functional seems to turn the phase transition into second order [4(|, in accordance with universality arguments. 

For phcnomcnological applications a parametrization of the model is needed, which is usually done at zero temper- 
ature. For this reason it is interesting to know which values of the parameters m 2 and A* at T* are of interest. In the 
region of the parameter space where there is a phase transition T c increases steeply with A*. Therefore, in order for 
the phase transition to occur at a reasonable value of temperature compared to T* one needs to restrict the possible 
values of the parameters in the to 2 — A* space. The accessible values of <j)(T = 0) and M(T = 0) when the restriction 
T c /T+ < 0.3 is imposed can be seen Fig. (JSJ. This means that the region of the parameter space which is interesting 
in practice is a relatively narrow one, close to the curve given in Eq. (fT2l . Setting T+ = 1 GcV, these values are in 
the range of the pion decay constant, pion masses, and the transition temperature in the 0(4) linear sigma model. 
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A*(m*) = 2m 2 



> Q 2 




FIG. 5: Dependence of the solution of the field equation at T = on the parameters in the region of the parameter space 
where T c /T* < 0.3 (T c increases with A*.) The labels near the solid curves refers to the solution M of the gap equation at 
T = 0. 



11 For practical purposes this curve can be approximated using the high temperature expansion: A*(m*) = 2m^/(T^/12 — m*T* /(87r)). 
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V. CONCLUSIONS AND DISCUSSION 

We have studied the phase transition of a real scalar (p 4 theory in four dimensions at lowest order in the 2PI 
formalism. Due to the existence of the Landau pole, this model needs to be considered in the presence of a cut-off 
lying below the scale of the Landau pole. For values of the parameters such that the Landau pole is much larger 
than all other scales, renormalization ensure insensitivity of the results with respect to the cut-off scale already below 
the Landau scale. We have illustrated these questions at the level of the gap equation by investigating under which 
conditions the cut-off can be continued above the Landau pole to infinity without noticing any sensible change in 
the solution of the gap equation. We have also studied the effective potential and proven analytically that if there 
is a phase transition, it is weakly first order, confirming earlier studies based on a numerical evaluation or a high 
temperature expansion of the effective potential. 

Let us end by emphasizing that even though we imposed the renormalization and consistency conditions in the 
symmetric phase, we could have imposed them in the broken phase as well. Because </> = is not always accessible 
we need to impose these conditions at <j) = v ^ (where v is arbitrary for the moment). The bare mass m 2 , and the 
bare coupling Ao are used to impose the renormalization conditions 

M|=„, t=o = ™? > and % =v ^ t=o = A > . (73) 

Unlike what happened at = in the symmetric phase, the curvature M^_ v T=0 is not equal to M 2 =v T=Q . This 
truncation artifact can be overcome by using the bare coupling A4 to impose the consistency condition: 

Mj, =v T=0 = M^ =v T=0 . (74) 

Notice then that the consistency condition: 

V ( f ) — V: t=o = y<i,=v, t=o j (75) 
fixes the value of v in terms of m? and A. An explicit calculation leads to: 12 

^ = f_ = 6 

~ m 2 A-3A 2 /(32tt 2 ) ' [ > 

One can then check that v cannot be interpreted as the value of </> (minimum of the potential) at T = which, given 
the renormalization and consistency conditions, is such that 4>x =0 /m 2 = <V^- Reversely, if one chooses v as the value 
of (f> at T — 0, which is the usual choice, then v 2 /m 2 = 3/A and one cannot impose the consistency condition (|75|) . 
In this later case, we have checked using similar arguments as above, that the transition cannot be second order. 

The impossibility to impose simultaneously the consistency conditions (|74[) and (|75[) at the minimum of the po- 
tential needs to be regarded as a truncation artifact. It can be cured by exploiting the fact that in higher orders 
of approximation of the 2PI effective action, two bare masses m 2 , and m 2 . are usually needed to renormalize differ- 
ent quadratic divergences in M 2 and M 2 , see for instance It is true that in the Hartree approximation these 
quadratic divergences are equal, but we could still allow for different finite parts to m 2 , and m 2 . and use them to 
impose both consistency conditions at the minimum of the effective potential. Notice that this approach and the one 
involving only one bare mass are not completely equivalent. They both agree at leading order in A but differ beyond 
this order: in the first approach, some of the consistency conditions, which reflect an exact property of the theory are 
sacrificed, whereas in the second approach, what is sacrificed is the fact that in the exact theory there is only one bare 
mass. The second approach can be viewed as a renormalization improvement over the first one, for exact properties 
of the theory are reflected at the level of the renormalized quantities not at the level of the bare ones. It is then 
interesting to wonder how such an improvement affects the order of the transition. Interestingly, before exploiting 
the consistency condition (|75|) a second order phase transition is possible for parameters such that v 2 /m 2 < 1/A. 
However, the consistency condition (|75[) leads once again to a relation between v 2 /m 2 and A: 

A 4 ( 3A \ i- 6 3A 3 v 4 v 2 „ 

1 )-K + 7^-T- A -o +3- (77) 



32tt 2 V 32tt 2 J m 6 16tt 2 m 4 



This shows in particular that A < 327r 2 /3 ~ 105. 
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For positive A, this equation has a positive solution v 2 /m 2 only for A < 327r 2 /3. Moreover, it is possible to show that 
v 2 /m 2 > 3/A and thus the region v 2 /m 2 < 1/A cannot be accessed. For parameters such that v 2 /m 2 > 3/A, it is 
again possible to convince oneself that the transition cannot be of second order. Then, our conclusions on the nature 
of the phase transition in the Hartree approximation remain unchanged even within this extended renormalization 
scheme. The impossibility of a second order phase transition is an artifact of the Hartree approximation. Therefore, 
it would be of interest to see to what extend one could apply the analysis and methods of the current work to more 
complicated approximations of the 2PI functional, in which there exist in the literature numerical evidence for the 
occurrence of a second order phase transition [40( . Investigations in this direction are under way. 

Acknowledgments 

We would like to thank J. -P. Blaizot, G. Fejos and J. Serreau for useful discussions on related topics. 

Appendix A: Useful identities 

With a three dimensional and rotation invariant regularization, the tadpole integral can be written as 

T - f l + 2n e f T * - f 6±rir 

G= - — - / G+ (Al) 



Q<A Jq<A 2e g JQ<A J q<A a q 

with e q = (q 2 + M 2 ) 1 / 2 , n s = l/(e l3e — 1) and £*n e = n e — n*, with n* the thermal factor at temperature T+. The 
infinitesimal version of Eq. (|Alj) is also useful and is obtained after taking a derivative with respect to T in Eq. (jAlj) 
followed by the limit T* — > T. We obtain 

r \ r T t\n /a r T \ 

(A2) 




with: 




(A3) 



and where the dot in h £ represents a derivative with respect to the explicit thermal dependence, so that n e = 
(e/T 2 )n e (l + n e ) > 0. Similar formulae can be obtained for the bubble integral, for instance: 

T & = f i + 2K- £ XJ = p &+ f (A4) 

?<A Jq<A Jq<a J q <A 2e g 



where the prime in n' E denotes a derivative with respect to energy, so that n' e = — (l/T)n e (l+n e ) < 0. The infinitesimal 
form of Eq. (|A4[) reads 

a ( r T \ r T Pini /a r T \ 

(A5) 




with: 




(A6) 



Appendix B: Solutions of the renormalized gap equation 

In order to discuss the solutions of the gap equation = <7a(M 2 ), we study the zeros of 

[M 2 )^-M 2 + ml + ^ ct> 2 + f T G- G* + (M 2 -ml) P G\ 

JQ<A JQ<A JQ<A 



(Bl) 
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A* 
2 



Q<A 



with first and second derivatives given by 

9a(M 2 ) = 

It is useful to notice that g\{M 2 ) and g' A {M 2 ) can be written as 



G 



G: 



Q<A 



A, / G 3 

Jq<a 



9a{M 2 ) 



g' A (M 2 ) = -A* 



(M 2 
1 



A* 
2 

G 2 



G 



Q<A 



Q<A 



Q<A 



(B2) 
(B3) 

(B4) 
(B5) 



We shall also use the explicitly renormalized form of (?a(M 2 ) given in Eq. (|2~Tj) . As explained at the beginning of 
Section [HI D( we perform our analysis in the regime of interest, that is where A and A p are much larger than all the 
other scales T*, m+, </>, and T. Moreover, we are only interested in positive solutions M 2 of the gap equation. 



Solution for A < A B 



The function g' A (M 2 ) is strictly positive, see Eq. (|B3[) . It follows that g' A (M 2 ) increases strictly from g' A (Q) = — oo 
(the bubble integral in Eq. (|B5jl diverges positively as M 2 — > 0) to g' A {oo) = — A*/Aq (the bubble integral in Eq. (IB5I) 



is suppressed as M 2 3> A 2 , T 2 ). Because A < A p , it follows from Eq. (|22|) and the definition of A p given in Eq. (|23|) 
that Ao > and therefore g^{oo) < 0. In consequence g A (M 2 ) < and the function g\(M 2 ) decreases strictly from 
«7 A (0) = A.(0 2 - C 2 (T, A))/2 where 



c 2 (T,A) 



2 mi 



Q<A 



1 

Q 2 



G+ 



Q<A 



T, 



G: 



(B6) 



<A 



to 9a(oo) = —oo (the tadpole integral in Eq. (|B4[) is suppressed as M 2 3> A 2 , T 2 ). Then, the existence of a solution 
to the gap equation depends on the sign of <M.(0) or, in other words, on the value of 4> 2 as compared to (f> 2 (T,A): if 
4> 2 > 4> 2 (T, A), the gap equation admits a unique and strictly positive solution; if <f> 2 = 4> 2 (T, A), the unique solution 
to the gap equation is equal to zero; if (j) 2 < 4> 2 (T, A), the gap equation has no solution. 

Notice that, so far, (j> 2 (T, A) does not need to be positive. In practice however, we shall consider only positive values 
of <fr and it is useful to determine when </> 2 (T, A) is positive as well. Using the explicitly convergent form of <?a(M 2 ) 
given in Eq. (|21[) . it is convenient to write 



0c (T, A) 



with 



= m 2 



A* 



2 a 

A, 



9 <A 9 



?<A 9 



Q<A 



Q 2 



(B7) 



(B8) 



If the parameters A, T*, m*, and A* are such that C* > 0, it is clear that 4> 2 (T, A) is negative for any value of T and 
therefore, it is of no relevance in practice because the gap equation admits a solution for any value of </> 2 > and T. 
On the contrary, when C* < 0, then 2 (T, A) becomes strictly positive for T <T C with the "critical" temperature T c 
defined by 



<z<A Q 



2C* 
A. 



(B9) 



In this case, 4> 2 {T,K) becomes relevant because the gap equation admits a (positive) solution only for cf> 2 > <fi 2 (T,A). 
Moreover, in the regime of interest, one can neglect the dependence of C* and of the thermal integrals with respect 
to A. The critical value cj) 2 (T, A) is then independent of A and one has (j) 2 c {T) = (T 2 - T 2 )/12 with T 2 = -24C*/A*. 
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2. Solution for A = A p 

From the definition of A p in Eq. (f23|) . this case corresponds formally to 1/Ao = 0. The function g' A (M 2 ) increases 
strictly from g' A (0) = — oo to g' A (oo) = — A^/A = 0. It follows that g' A (M 2 ) < and the function g A (M 2 ) decreases 
strictly from g\(0) = A*(</) 2 — 4> 2 {T))/2 to g A (oo) = A*(0 2 — 0p)/2 (the tadpole integral in Eq. (|B4[) is suppressed as 
M 2 > A 2 ), with 

<f>l = [ * G* > . (BIO) 

JQ<A P 

The existence of a solution to the gap equation now depends on the signs of <?a(0) and g\(oo) or, in other words, 
on the value of <j> 2 as compared to <j> 2 (T) an <f> 2 . Notice that, because <?a(0) > g A (oo), we have <fi 2 > 4> 2 {T). We can 
thus distinguish three different regimes: if (f> 2 > cf> 2 , the gap equation has no solution; if <fi 2 > <fi 2 > <j> 2 (T), the gap 
equation admits a unique solution, which vanishes for <fi 2 = <fi 2 (T); if <fi 2 < (j) 2 (T), the gap equation has again no 
solution. Notice that, in the regime of interest, 2 ~ A 2 /(87r 2 ) is a large scale. 

It is important to stress that the behavior of the solution M 2 of the gap equation as A — >• A" is rather different 
depending on the value of (j> 2 . To see this, note that for a given /j,, the limit of <7a(m 2 ) as A — > Ap is A*(0 2 — <fi 2 + 
Iq<a Gfj)/2, sec Eq. (|B4[) . Clearly, if (f> 2 > <p 2 , this limit is strictly positive for any /i, which means that for any 

[/,, there exists a value of A < A p above which g^u 2 ) > and thus M 2 > fi 2 . This shows that the solution of the 
gap equation diverges as A — >• A" when <\> 2 > <fi 2 . If <f> 2 < <fi 2 , nothing of this kind happens and the solution can be 
extended beyond the Landau scale. We now study what happens beyond this scale. 

3. Solutions for A > A p 

This case corresponds to Ao < 0. The function g' A (M 2 ) increases strictly from g' A (0) = — oo to g' A {oo) = — \+/\q > 0. 
It follows that g' A (M 2 ) vanishes for some M 2 (T, A) (in view of Eq. (|B5j) this mass does not depend on <f> 2 ). Then, 
the function g A (M 2 ) has a minimum at M 2 = M 2 (T,A): it decreases strictly from g\(0) = A*(</> 2 — 2 (T))/2 to 
g^{M 2 {T, A)) = \*(<j) 2 — 2 (T, A))/2 and then increases towards 5a(oo) = oo (the tadpole integral in Eq. (|B4|) 
is suppressed as M 2 ^ A 2 ). The existence of solutions to the gap equation depends now on the signs of ,9a (0) 
and gA(M 2 (T, A)) or, in other words, on the value of 4> 2 as compared to (f> 2 (T) and tfi 2 (T,A). Notice that because 
<7a(0) > g\(M 2 (T, A)), we have 4> 2 (T,A) > 4> 2 (T). We can thus distinguish four different regimes concerning the 
number of solutions to the gap equation: if <f> 2 > <pp(T, A), there is no solution; if <fi 2 = 4>p{T, A), there is one solution, 
equal to M 2 (T, A); if </> 2 (T, A) > 4> 2 > 4> 2 {T), there are two solutions, one smaller than M 2 (T, A), the other larger 
than M 2 (T, A); if 4> 2 < 2 (T), there is again one solution, larger than M 2 (T, A). 

It is useful to study the dependence of M 2 (T, A) and </> 2 (T, A) on A. Let us consider M 2 (T, A) first. It is defined 
by g' A (M 2 (T, A)) = 0. Using Eq. (|B5|, it follows already that M 2 (T, A) ->• oo as A -> A+, because -1/A -4- 0+. Let 
us now show that, in the regime of interest, M 2 (T, A) > m 2 or, in other words that g' A {m 2 ) < 0. From Eq. (|B2[) . we 
write 



9A( m l) = 


-1 - 


A, ' 
2 


f gi - r gi 

JQ<A JQ<A 


< 


-1 - 


A* ' 
2 


r T=0 r T t 

/ Gl- Gl 

J Q<A JQ<A 


< 


-1 + 




r nt* — sin*'* 




2e* 3 ' 



(Bll) 



where we have used Eq. (|A4[) and the fact that the bubble integral increases with T and A. Then <? A (m 2 ) < if we 
require that 
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In terms of the Landau pole this condition reads 



16tt 2 



Arcsinh 



A, 



A, 



A 2 + ml 



> 



-q' 1 e{ 



g>A p 



4e* 



(B13) 



which is clearly fulfilled in the regime of interest. Finally, in order to study the A-dependence of M 2 (T, A), let us take 
a derivative with respect to A in the defining equation = g' A (M 2 (T, A)). We arrive at 



1 



dg' A 



dA gl{MD dA 
Using Eqs. (|B2[) and (|A4[) . this reads more explicitly: 



dMl 
dA 



A. 



A- 



2^(M p 2)2^ 



1 + 2(rc,p 



e p A n' eP ) 



4 ,r 



1 + 2(n* x e\n% 
4 £ * 3 



(B14) 



(B15) 



Now, because T <T+ and Mp(T, A) > m 2 , the right-hand-side of this last equation is strictly negative and Mp(T, A) 
decreases strictly with A. Because it is bounded from below, it has a limiting value M 2 (T, oo) which one could 
determine from Eq. (|B2[) . In the regime of interest, one obtains Mp(T, oo) ~ 4A 2 /e 2 . 

Consider finally 4>Z,(T, A) which is defined as the particular value of 4> 2 such that g\(M 2 (T, A)) is equal to zero. 
Making the field dependence of g A (M 2 ) explicit we can then write 

= . 9A (Af p 2 (T,A),^(T,A)). 
Taking a field derivative with respect to A and using the fact that g' A (M 2 (T, A)) =0, we obtain: 







dgAd&l dg A 



dA dA 

A straightforward calculation using the explicitly convergent form of g A (M 2 ) given in Eq. (|2 1 1) leads to: 

A 2 



9% 
dA 



2- 



,2\2 



G; G p + thermal . 



(B16) 
(B17) 

(B18) 



where G p = 1/(Q 2 + M 2 ), the propagators are evaluated for q 2 = A 2 and q\ denotes the continuous Euclidean 
frequency at T = 0. In the regime of interest M 2 (T, A) > M 2 (T, oo) S> to 2 - One can thus neglect the thermal 
contribution in the previous formula. It follows that 4&(T, A) decreases strictly with A. For A — > A+, the minimum 
M 2 (T,A) of g A (M 2 ) is sent to oo. Thus <?a(A? 2 ) = A*(</> 2 - c/> 2 (T,A))/2 -> g Ap (oo) = A*(> 2 - c/> 2 )/2 and thus 
</> 2 (A+) = 2 Moreover, using Eqs. dHTJ) and dB2j) with g'^M 2 (oo)) = 0, we obtain 



^p(oo) 



5*n e 



-± + (M^(oo) - miY G 2 (oo)G* 



(B19) 



Q 



which is strictly positive for T <T±. In the regime of interest, only the zero temperature contribution of the second 
term dominates and we find 



#(«>) 



AJp(oo) 
8tt 2 



\ 2 
p 



(Q 2 + l) 2 (Q 2 + m 2 /M 2 (oo)) 47r 2 e 2 



Notice finally that : 



dM 2 



A 



(B20) 



(B21) 



dej? 2 g' A (M 2 ) ■ 

This shows that, as one increases (f> 2 , one of the solutions increases while the other decreases. They both become 
equal to M 2 (T, A) when (f> 2 = </> 2 (T, A). This applies in particular to the case A = oo. 
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